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Abstract 

In this paper, we study the loss of coherence of a wave propagating according to the Schrodinger 
equation with a time-dependent random potential. The random potential is assumed to have 
slowly decaying correlations. The main tool to analyze the decoherence phenomena is a prop- 
^N) erly rescaled Wigner transform of the solution of the random Schrodinger equation. We exhibit 

^-H anomalous wave decoherence effects at different propagation scales. 
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Oh 1 Introduction. 
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a 



The Schrodinger equation with a time-dependent random potential has attracted lots of attention 
because of its large domains of apphcations |T5J [T21 EH HSl 137] • It is widely used for instance in 
wave propagation under the paraxial or parabolic approximation 131 HI El |7] . This field of research 
was recently stimulated [71 UHl 1221 122] EHl ES] by data collections in wave propagation experiments 
showing that the medium of propagation presented some long-range effects [nj ES]. Most of the 
theoretical studies regarding wave propagation in long-range random media hold in one dimensional 
^ propagation media, which are very convenient for mathematical studies but not relevant in many 

applications. 

In this paper, we consider the following random Schrodinger equation 

^ 1 

idt(t>+ -A^^- ^V{t,^)(j>^0, t>OandxeM'*, 

^ 0(O,x) = 0o(x), 

". I with a random potential V{t, x), which is a spatially and temporally homogeneous mean-zero random 

J>] field. Here, t > represents the temporal variable, x S M'* the spatial variable with d> 1, and e <C 1 

is a small parameter which represents the relative strength of the random fluctuations. A classical 
tool to study the loss of coherence of the wave field is the Wigner transform defined by 



X 



which is somehow the Fourier transform of the correlation function in space of (j> at point x. The 
Wigner transform measures the degree of correlation of the wave field </> at two different points (x— y/2 
and X -|- y/2). The loss of coherence of the wave field ip (or equivalently the wave decoherence) 
corresponds to the evolution of the degree of correlation in space of the wave field (j) and is captured 
by the wave number k through the evolution in time of the momentum of the Wigner transform W. 
Let us remark that if y = in ([T|), we have M^(i,x,k) — IVo(x — tk,k), where Wq is the Wigner 
transform of the initial data 0o- In this case the momentum of W is preserved during the propagation. 
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that is there is no variation of the momentum with respect to time meaning that there is no loss of 
decoherence. 

We refer to [HJ [57] for the basic properties of the Wigner transform. In our problem the amplitude 
of the random perturbations are small, so that to observe significant cumulative stochastic effects we 
have to look at the wave field </> over long times and large propagation distances. Consequently, we 
consider the rescaled field 

0,(i,x) = 0(^,ii) (2) 
satisfying the scaled random Schrodinger equation 



^2s 



ie'dtde + VAx^e - \^v( — , = 0, t > and X G R'^, 

2 \e'^ / 

0£(O,X) = (?!)o,e(x), 

and where s 6 (0, 1] is the propagation scale parameter. If the random potential V has rapidly 
decaying correlations, it has been shown [7] that the wave field (/)<: does not exhibit any significant 
random perturbations before the propagation scale (s = 1). More precisely, the Fourier transform 
in space of 0^ (defined by ([2|) properly scaled converges point- wise to a stochastic complex Gaussian 
limit with a one-time statistic of an Ornstein-Uhlenbeck process. Wave decoherence phenomena have 
been studied in this context in many papers [3J|ll[I21[I31[ni[ISl[I71[231llHlE7] thanks to the following 
Wigner transform for s = 1 



W^,(i,x,k) = — . / dye^''-y0eU,x-e^UeU,x 



{2ttY J V 2/^'V ' 2 



e- 



(27r)° 



(3) 



The Wigner transform ([s]) analyzes the loss of coherence of 4>e on a spatial correlation scale of order e 
at the macroscopic propagation scale e~ ^ , and the loss of coherence of on a spatial correlation scale 
of order 1 at the microscopic propagation scale. Let us note that 



\ct>,{t,x)\'' = J dW,(t,x,k), 



so that we may think of as a wave number resolved energy density. However, the terminology 
density is not well appropriate since is not always positive but it becomes positive in the limit 
e^O [27]. 

In several context involving rapidly decorrelating random potentials, it has been shown that the 
expectation of the Wigner transform E[VF(:(<:, x, k)] converges as e goes to to the solution W of the 
radiative transport equation 

dtW + k-V^W = y"dp(T(p,k)(T4^(t,x,p)-T4^(t,x,k)), (4) 

where the transfer coefficient cr(p, k) depends on the power spectrum of the two-point correlation 
function of the random potential V. Moreover, in some cases [2J[11[^[TB], it has been shown that the 
limit W is often self-averaging, that is We converges in probability to the deterministic limit W for 
the weak topology on L^(K^'^). In other words, the wave decoherence mechanism described by the 
radiative transfer equation Q does not depend on the particular realization of the random medium. 

In this paper, we investigate the loss of coherence of a wave propagating according to the Schrodinger 
equation ([ij involving a random potential V with slowly decaying correlations defined in Section 2.1 



In this context, the behavior of the scaled field (f>t defined by Q evolves on different propagation 
scales [7] [22]. In [T] Theorem 1.2] the authors study (pe itself on the propagation scales e~*, 
with s = 1/(2k) and k > 1/2, and show that the Fourier transform of 0e properly scaled converges 
point-wise in distribution to a complex exponential of a fractional Brownian motion with Hurst index 
K, where k depends on the statistic of the random potential (see Section |3| . This result means that 
e^'^, with s = 1/(2k) < 1, is the first propagation scale on which the random perturbations become 
significant, and induces a random phase modulation on the propagating wave. In |22l Theorem 2.2] 
the author studies the loss of coherence of (j)^ on the propagation scale e^^ (s — 1), and shows that 
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the Wigner transform (|3| of 0^ converges in probability for the weak topology on to the 

unique solution of a deterministic radiative transfer equation similar to Q. In other words, the wave 
decoherence happening on this propagation scale (e~^) does not depend on the particular realization 
of the random potential. However, even if the radiative transfer equation has similar structure in 
both cases (for rapidly and slowly decaying correlations), the presence of a potential with long-range 
correlations have a striking effect. In contrast with the rapidly decorrelating case studied in the 
scattering coefficient S(k) = J (ipcr(k,p) = +00 is not defined anymore. The radiative transfer equa- 
tion is however still well defined because of the difference W{t,x,p) — W{t,x,'k) which balances the 
singularity introduced by the long-range correlation assumption. Moreover, these long-range correla- 
tions imply an instantaneous regularizing effect of the radiative transfer equation |22| Theorem 3.1]. 
Consequently, these results show a qualitative and thorough difference between the rapidly and slowly 
deccorelating cases. In fact, in contrast with the rapidly decorrelating case, for which the phase and 
the phase space density evolve on the same propagation scale [7j, the phase of 0e and its phase 
space energy now evolve on different propagation scales. 

The main goal of this paper is to study the loss of coherence of a wave propagating according to 
the Schrodinger equation (jlj involving a random potential V with slowly decaying correlations for 
propagation scale parameters s € (1/(2k),1). In fact, for s > 1/(2k) the random phase modulation 
obtained for s = 1/(2k) implies that the wave field cf)^ produces very fast phase modulation, so that 
for sufficiently long times and large propagation distances e^* (with s > 1/(2k)) the wave coherence 
should be broken. Let us note that for a given propagation scale parameter s, the wave decoherence 
can be too small to be observed. For instance in [22 , the author shows using the Wigner transform 
([3]) that there is no significant wave decoherence on the wave field on spatial scales of order e*, 
before s = 1 . As we will see in Section |4] for s < 1 wave decoherence takes place first on large spatial 
scales and then propagates to the smaller ones as the propagation scale parameter s increases. The 
larger the propagation scale parameter s is the smaller the spatial scale is to observe wave decoherence 
(see Figure [2| . To exhibit wave decoherence for s < 1 , we need to consider a properly scaled Wigner 



transform of the field 0^ (see Section 2.3). Depending on the propagation scale parameter s, we show 
that this scaled Wigner transform converges in probability, for the weak topology on L^(M^'') to the 
unique solution of a fractional diffusion equation. This momentum diffusion equation describes the 
wave decoherence mechanism, and shows that it does not depend on the particular realization of 
the random potential. The anomalous momentum diffusions obtained for s £ (1/(2k), 1] allow us to 
exhibit a damping coefficient, describing the decoherence rate, obeying a power law with exponent in 
(0,1). 

The organization of this paper is as follows: In Section |2] we present the random Schrodinger 
equation that will be studied in this paper; then we present the construction of the random potential 
and we introduce the long-range correlation assumption on the potential V used throughout this 
paper; finally, we introduce the rescale Wigner transform which is the main tool to describe the 
wave decoherence mechanisms. The results stated in Section [3] and Section |5] have been respectively 
shown in [7] and , but we recall these results to provide a self-contained presentation of the wave 
decoherence phenomena. In Section [3] we present the behavior of the field 0e for the propagation 
scale parameter s = 1/(2k). In Section |4] we state the main result of this paper. We present the 
asymptotic behavior in long-range random media of a properly scaled Wigner transform over the 
intermediate range of propagation scale parameter s € (1/(2k),1). In Section |5] we describe the 
asymptotic evolution in long-range random media of the phase space energy density of the solution 
of the random Schrodinger equation for s = 1. Finally, in Section [6] we recall the proof of Proposition 
\2.1\ and Section [7] is devoted to the proofs of Theorem |4.1| and Theorem |4.2[ 
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2 The Random Schrodinger equation 

This section introduces in a first time the random Schrodinger equation studied in this paper, then in 
a second time it introduces the random potential with long-range correlation properties. Finally, we 
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introduce the Wigner transform which is the main tool in this paper to study the loss of coherence of 
solutions of the Schrodinger equation. 

We consider the dimensionless form of the Schrodinger equation on M'' with a time-dependent 
random potential: 

z9t</)+^A,,^~e^F(^,x)0-O, (5) 

with 76 [0, 1). 7 is a parameter characterizing the correlation length in time. If 7 = the correlation 
lengths in space and time are of the same order, but if 7 g (0,1) the correlation length in time is 
small compared to the correlation length in space. In ([s]) the strength of the random perturbations 
are small, so that to observe significant cumulative stochastic effects we have to look at the wave field 
(j) over long times and large propagation distances. Consequently, throughout this paper we consider 
the following rescaled wave field : 

0,(i,x) = J), with se(0,l]. (6) 

Moreover, throughout this paper the parameter s £ (0, 1] represents the propagation scale pa- 
rameter, and the scaled wave field 0e satisfies the scaled Schrodinger equation 



+ ^Ax0.- 6^^(^,^)0,-0 with (/.,(O,x) = 0o,.(x). (7) 

Here Ax is the Laplacian on M'* given by A = '^xj ■ (^(^j x), i > 0, x G M"*) is a random potential 

given by a stationary zero-mean continuous random process in space and time, and whose properties 
are described in Section 2.1 The initial datum 0o.e(x) = <^o^e(x, is a random function with respect 
to a probability space (5*, 5, /Lt(dC)), and independent to the random potential V. This randomness 
on the initial data is called mixture of states. This terminology comes from the quantum mechanics, 
and the reason for introducing this additional randomness will be explained more precisely in Section 

[231 

2.1 Random potential 

This section is devoted to the introduction of the random potential V considered in this paper, and 
is also a short remainder about some properties of Gaussian random fields that we use in the proof 
of Theorem |4.1| and Theorem |4.2[ All the properties of the random field V exposed in this section 
result from the standard properties of Gaussian random fields presented in [T] for instance. 

Before introducing the random potential we need some notations. Let Rq be a positive function 
such that i?o G L^(M'^), _Ro(— p) = ^o(p), decaying rapidly at infinity and having a singularity at 
p = 0. Let us consider the following Hilbert space : 

'Hm = i 'Z' such that <p(p) = </5(— p) and / TO(dp)|(^(p)|^ < -l-oo >, 

equipped with the inner product 



{v^'^)-H "i(c^P) ¥'(P)^(P) V((p,'i/') e H„i X "H™, 



and where m{dp) — Ro{p)dp is a real finite measure. Let (wn)n>o be an orthonormal basis of Jir 
and let (Vn)n>o be a family of real-valued stationary zero-mean Gaussian processes such that 

HVp{t)%is)]= [ TO(dp)e-8(P)l*-"lup(p)uq(p), Vs,i > and Vp,g > 0. 
Now, let us consider the following real-valued linear functional on Hm defined by 

" n>0 

where 'H^ stands for the dual space of Hm , which is well defined since 



" p,q>o 
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As a result, (^^(i))t>o is a real- valued stationary zero-mean Gaussian process on T-L'^ such that 

for all t, s > and for all ((/j, ^) e 'Hm x Hm, which corresponds to a covariance function given by 

nV{t, dpi)F(s, dp2)] = (2^)'^i?(i - s, pi)<5(pi P2), 

where the spatial power spectrum is given by 

i?(t,p)-e-8(p)l*li?o(p). (8) 

Here, the nonnegative function g, such that g(p) — fl(~p), is the spectral gap. Particular assump- 
tions involving the spectral gap q will be introduced at the end of this section to ensure long-range 
correlation properties on the potential V in ([t]) which is defined as follows : 

V{t, x) = I V{t, dp)e'P-^, Vt > and Vx G , (9) 

(27r)'' J 

so that (y(i,x),t > 0,x e W'-) is a real- valued stationary zero-mean real Gaussian process with 



E[\/(t,x)F(,s,y)] =i?(i-.s,x-y) = ^ j dpR{t ~ s,p)e'P<^-y^ 

(^^^ J dujdpR{u;,p)e 



^ ' (10) 



for all t, s > and for all (x, y) G M~ , and where space and time power spectrum is given by 

^(->P)=^#4rT- (11) 
Moreover, according to [1, Theorem 1.4.1 pp. 20] and the fact that 

E[{V{h,x) - V{t2,y))Y^^ < C (^1 dpi?o(p)) (|ii-t2| + |x-y|), V(ii,t2,x,y) € [0,T]^ x 

the random potential V defined by ([9| is continuous and bounded with probability one on each 
compact subset K of M+ x M''. Finally, according to the shape of the spatial power spectru m ([s] ), we 
have the following proposition which will be used in the proof of Theorem 4.1 and Theorem |4 . 2 1 based 
on the perturbed-test-function method. 

Proposition 2.1 Let 

Tt = aiVis,-),s<t) (12) 
be the a -algebra generated by (V{s, ■), s <t). We have 

E[V{t + h,-)\^t] = e-8(P)''i>(i, •) (13) 

and for all {(p, tp) G Hm x Hm 



E 



(Vit + h),^)^,^^^JV{t + h),4')^,^^^^ -E[{V{t + h),^)^.^^.^^^ \Tt]E[{V{t + h),i;)^,^^^^ \Ft] 

= j dp ^(p)V(-p)i?o(p) (1 - e-^sCp)/.) . 

(14) 

The proof of Proposition |2 . 1 1 is standard but we still prove it in Section [6] just to show how to obtain 
these results in a weak formulation. 
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2.2 Slowly decorrelating assumption 

In this paper we are interested in the Schrodinger equation with a random potential with long-range 
correlations. Let us introduce some additional assumptions on the spectral gap q of the spatial power 
spectrum ([s]) in order to give slowly decaying correlation properties to the random potential V defined 
by ^. 

Let us note that for all fixed t > 0, the random field V{t,-) has spatial slowly decaying correlations. 
In fact, if we freeze the temporal variable, the autocorrelation function of the random potential V(t, ■) 
is given by 

i?(i,x) =E[^^(t,x + y)l/(i,y)] = ^ J dpi?o(p)e'P-'' 

where i?o(p) is assumed to have a singularity in 0, so that R{t, ■) ^ L^{M.'^). As a result, (V{t))t>o 
models a family of random fields on M.'^ with spatial long-range correlations which evolves with respect 
to time. However, since ([s]) is a time evolution problem, we have to take care of the evolution of 
the random potential V with respect to the temporal variable. In fact, if V has rapidly decaying 
correlation in time, {V{ti),V{t2)) has now rapidly decaying spatial correlations, and the evolution 
problem ([S]) behaves like in the mixing case addressed in [B]. As a result, even if at each fixed time 
the spatial correlations are slowly decaying, the resulting time evolution problem behaves as if the 
random potential has rapidly decaying correlations. Consequently, we have to introduce a long-range 
correlation assumption with respect to the temporal variable. 

For the sake of simplicity, throughout this paper we assume that 

fl(p)^HpP^ and flo(p)= |pff_i) ^ (15) 

where a is a positive continuous function decaying rapidly at infinity and such that a(0) > 0. Moreover, 
we assume that /3 S (0,1/2], a G (1/2,1), and a + (3 > 1, so that 

^ -^o(p) , 

dp — ^— = +00, (16) 
0(P) 

since 

^o(p) a(0) 



0(P) P-o|pr« ^ = 2(a + /?-l)E(0,l). (17) 

A similar configuration has already been considered in [7. to study the propagation of the wave field 
(j), defined by ^ in a random media with long-range correlations. As a result, these assumptions 
permit to model a random field V{t,ii) with spatial long-range correlations for each fixed time t>0, 
and with slowly decaying correlations in time because of the following relation : 

V(s,x,y) e M+ X M^d r dt\E\V{t + s,^ + y)V{s,y)]\=+oo^ f dp^^:^ = +oo. (18) 
Jo J 0(P) 



If fact, if (16 1 holds, we have 



lim / dt\E\V{t + s,ii + y)V{.s,y)]\> lim / dp^^^^(l - e-^'P)-^) - / dp^^^^je^P'" - 1| 
A^+^Jo ' ^ ' '-"'Jl - A^+oo7 ^ fl(p) ' J 0(p) ' 



and the converse implication is obvious by taking x 0. Consequently, throughout this paper we say 
that the family (V{t))t>o of random fields with spatial long-range correlations has slowly decaying 
correlations in time if (|16[) holds, and rapidly decaying correlations in time otherwise. 



2.3 Wigner transform 



In this paper we study wave decoherence phenomena, using the following Wigner transform ( 19 1 of the 
wave field ^ satisfying the Schrodinger equation ([T]), happening on different propagation scales s. In 
this paper we consider the Wigner transform of the field , averaged with respect to the randomness 
of the initial data, defined by: 



W^,(t,x,k) = -^ / dy/i(dC)e^'''^0,(t,x-e^-^=^,C)0.(t,x + e^-^ey ^) 
(27r)''7Rdxs V 2 y V 2 J 



(19) 
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where Sf. £ [0, s] is the spatial correlation parameter, and {S,S, jj.{d()) is a probability space. 



Let us remark that the scaled Wigner transform (19 1 is a real-valued function. We discuss below 



the reason of introducing this probability space, and we refer to [27^ for the basic properties of 



the Wigner distribution. The scaled Wigner transform (19 1, which is somehow the Fourier transform 
of the correlation function in space of the wave field (j) (resp. (f>e) around x, measures the degree of 
correlation of the wave field. In other words, it captures the evolution of the degree of correlation of 
the wave field (j) over a spatial correlation scale of order e^^", on the microscopic propagation scale. 
Or in the same way, it captures the evolution of the degree of correlation of the rescaled wave field (f)^ 
over a spatial correlation scale of order e*^*'=, on the macroscopic propagation scale e~*. 

The loss of coherence of the wave field satisfying ^ (or equivalently the wave decoherence) 
corresponds to the evolution in time of the momentum of the Wigner transform W^. In other words, 
if the momentum of the Wigner transform is preserved during the propagation there is no wave 
decoherence, and in the opposite case the wave decoherence mechanism is described by the evolution 
of the momentum. 

In this paper we consider a rescaled version of the Wigner transform involving the two parameters 
s and Sc- The propagation scale parameter s has been introduced in Section [l] and characterizes the 
order of magnitude of the propagation times and the propagation distances on which we consider 
the wave field satisfying (|5]). The spatial correlation parameter Sc characterizes the spatial scale 
on which we study the evolution of the degree of correlation of the wave field (p. Let us note that 
the cases Sc < s study the local loss of coherence of the wave field, while the case s — Sc study 
the nonlocal loss of coherence of the wave field. We will see in Section |4] and Section [5] that for a 
given propagation scale parameter s the loss of coherence of the wave field (j) can be observed at a 



particular spatial correlation parameter Sc depending on s (see (28)). The larger the propagation 
scale parameter is the shorter the decoherence scale parameter is. In other words, depending on the 
propagation time and propagation distance we adjust the spacial correlation scale characterized by 
Sc to exhibit the loss of coherence (see Figure |2]). For instance, according to |22] no significant loss of 
coherence of the wave can be exhibited on the correlation scale = before s = 1. Let us note that 
for a rapidly deccorelating potential V , no wave decoherence effects can be observed except for the 
radiative transfer scaling Sc = and s = 1 |7j . The reason will be explain formally in Section |4] 

However, to observe decoherence effects of the field (pe on the spatial correlation scale of order 
e*-^<:, we need a proper initial condition (t>o^f_ in ([7|, which oscillates at the same scale (see Figure 
[ij. Moreover, a natural way to introduce randomness on the initial condition is as follow. Let 
S — M'* and /i(C) be a nonnegative rapidly decreasing function such that ||/i||Li(R'') = li s-nd so that 
(M'',S(M'^'), /i(c?C)) is a probability space. Throughout this paper we assume that the initial condition 
00,6 in ^ is given by 

0o,e(x) = 0o(x) exp(iC • x/e^-^^). (20) 

This initial condition represents a plane wave with initial propagation direction C, G R'^, oscillating 
on the spatial scale e''"'''^, and with amplitude or envelope (pQ. The initial direction ( of the wave 



is distributed according to fJ,{d(), so that the Wigner transform (19 1 is average according to the 
distribution of the initial direction of the wave. Let us note that the spatial frequency of the initial 
condition (^ g~('*^''<:)) is low compared to the one of the random medium (^ e^**) on the macroscopic 
scale e^^. In rapidly deccorelating random media such low spatial frequency sources do not interact 
with the random medium, but as we will in Section [4] this kind of initial conditions interact strongly 
with slowly decorrelating random media. This result can be useful in passive imaging of a target in 
slowly deccorelating random media [T^. 

The main reason to introduce this additional randomness through the initial data 4>o.e is to make 
possible the weak convergence of the initial Wigner transform W^c(O) in L^(M^''), which stands for the 
set of real-valued square-integrable functions equipped with the following inner product 

(A,/i)^ = / dxdk/(x,k)5(x,k), V(A,/i)eL2(ig2<i)^^2(i^2d)^ 

Consequently, we have 

VA e L^R^^), lim (14^,(0), A)^,(j^,,^ = (Wo, X) 

where 

W^,(0,x,k) = W^o,.(x,k) = / dyM(dC)e*''^-'^)-'^0o(x - e^-^=y/2)0o(xT^^=^^, (21) 
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-5 -1 -3 ^ 




Figure 1: Illustration of the initial condition (20 1, (a) and (6) represent low spatial frequency initial 
conditions compared to the spatial frequency of the random medium e"" on the macroscopic scale 
e~*. (a) represents the case Sc = s and has a spatial frequency of order 1 on the macroscopic scale 
(and a spatial frequency of order on the microscopic scale) . (6) represents the case Sc < s and has 
a spatial frequency of order e"'-""*'''' on the macroscopic scale (and a spatial frequency of order 
on the microscopic scale), (c) represents the case Sc = and s = 1, and has a spatial frequency of 
order on the macroscopic scale (and a spatial frequency of order 1 on the microscopic scale). 
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with 

W^o(x,k) = |0o(x)|V(k), ifsc<s, (22) 

and 

M/o(x, k) = -i-^ / dyfiidOe'^'^-^yyM^ - y/2)0o(x + y/2), if s = s,. (23) 

Consequently, thanks to the Banach-Steinhaus Theorem, Wo,^ is uniformly bounded in L^(M^'') with 
respect to e. We need such a convergence on the initial Wigner transform M^e(O) since we study the 
Wigner transform (19 1 in in i^(E^'') equipped with the weak topology. As it will be discussed 
in Section [4] it is not possible to expect a convergence result in L'^{MP'^) equipped with the strong 
topology (except for the case s = Sc)- 

The Wigner distribution ( 19 1 satisfies the following evolution equation 

dtW.it, X, k) + e^^k • VMt, X, k) = 

V(j^,dp) ..... p . . p (24) 



g(l-7)/2-s 



with initial conditions VFe(0,x,k) = M^o.e(x,k), and where Wo^e is defined by (211. The previous 
equation ( 24 1 can be recast in the weak sense as follows : 



{W,{t),X)^,^^,,^-{WM,^)mR2.-, = (W-.(w),e^=k.VxA + e(i-^)/2-^£,A(^.))^,(^,,) 

for all A e 5(M2rf), where ^(M^rf) is the Schwartz space and stands for the space of rapidly decaying 
functions. Here 



^^^(^'-'"^^ = (2^ 2^) 2^ 



(25) 



Let us assume that y = 0, so that the Wigner transform is given by W^{t, x, k) = H/o(x — e*^ik, k). 
Therefore, the momentum of Wc is preserved during the propagation, meaning there is no variation of 
the momentum with respect to time, that is there is no wave decoherence. The dispersion term k • Vx 
is of order e^" so that captures the wave dispersion only for Sc = 0. The transfer equation (24 1 
describes the loss of coherence of the field (f>^ through the random operator £eWe(i)(t, x, k). However, 
depending on the spatial scale of observation the loss of coherence of the wave may not be significant. 
In fact, according to no significant wave decoherence can be exhibited on the correlation scale 
Sc = before s — 1. 

Let us introduce some notations which are used in Section |4] and Section [S] Let 

Br = {Xe L^{R^'^),\\X\\L2(M^<i) <r} , with r = sup || W^o.e||L2(R2d) < +00, 

be the closed ball with radius r, and {gn,n > 1} be a dense subset of Br- We equip Br with the 
distance dg^ defined by 



+00 ^ 

rfe,(A,Ai) = X! ^ |(-^ ~ ^'5")l2(r2'') ' V(A,Ai) e 6,. X 6^, (26) 

so that (Brjdg^) is a compact metric space. Therefore, (W^c)c is a family of process with values 
in {Br,dB^), since || W(:(i)||i2('K2d') = || Wo,e||i2(j{2d'). The topology generated by the metric dg^ is 
equivalent to the weak topology on L^(M^'*) restricted to Br- 

The three following sections describe in a chronological order the effects produced by the random 
medium on the wave propagation. 



3 Phase Modulation Scaling s = 1/{2k,^) 

This section describes the first effects caused by the small random fluctuations of the medium on 
the wave propagation. The following theorem presents the asymptotic behavior of the phase of (/)e 
solution of ([7]). Theorem 3.1 has been shown [7^ in the case 7 = and Sc = 0, but nevertheless, its 



proof remains the same as the one of [7) Theorem 1.2]. We state this result in order to provide a 
complete and self-contained presentation of the wave propagation in long-range random media. 
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Theorem 3.1 Let us note 



and 



V 7 6 [0,1), 



l-7( 

and let us consider the process C.K^,e{t,'k) defined by 

C«„.(i,k) = -^^^0,(i,-j^)e^l'^l'*/(2.°""^), with s = 1/{2k^), ands,<s, (27) 

where 0e satisfies ([T]) with initial data (20). Under the long-range correlation assumption in time 
(16 1, where a + /3 > 1, the process CK^.e(i,k) converges in distribution to 

C(t,k) = Co(k)exp (^v/^(a,/3,k)B,„(^)), 

for each t > and k € M.'^, where (i?Ko(i))t>o ^'^ i standard fractional Brownian motion with Hurst 
index kq, and 

Co(k) = (f'aiC ~ k) if Sc = s, and Co(k) = (/'o(0)(J(C ~ k) otherwise. 

Moreover, 

i?(a,/3,k)=— -^AJ dp^^l dS{u)e\'^\P^-'^- if/3= j^0,ands, = 0, 

(27r)''Ko(2Ko - 1) Jo p^" ^ Jsd-i 2 



and 



i?(a,/3,k) = D(a,/3) = ^^^^^^^-^^ .^W.e, 



where, i^d is the surface area of the unit sphere in M.'^, and ei G S'' 



In Theorem 



3.1 



Co represents the initial direction of the wave, and exp(i-y/D(a, /?, k.)B^g (t)) 
represents the random phase modulation induces by the slowly decorrelating perturbations of the 
medium. The scaling k/e''"'*'^ in (27 1 corresponds to the order of the spacial frequency of the initial 
condition ( 20 1 . The order of the spacial frequency of the initial condition does not play a significant 
role in Theorem 13.11 but it will in Section H] and Section [H In fact in these sections the order of the 
spacial frequency of the initial condition (pO]) plays a significant role and is linked to the correlation 



scale parameter Sc of the Wigner transform (19 1 to exhibit the loss of coherence. The order of the 



spacial frequency is taken into account in Theorem |3.1| in order to give a complete picture of the wave 
propagation phenomena (see Figure [2]) . 

As a result, in long-range random media macroscopic effects may happen on the field 4>e at a 
shorter scale s = 1/(2k^) < 1 without induced loss of coherence of the field (p^. These effects are just 
a phase modulation given by a fractional Brownian motion with Hurst index kq, which depends on 
the statistic properties of the random potential V. However, let us note that the propagation scale 
parameter s — 1/(2k^) < 1 is "universal" in the sense that a random phase modulation appears on 
the wave whatever the order of the low-frequency initial condition characterized by Sc, that is for all 
Sc e [0,s = 1/(2k^)]. 

[7] is the first paper showing a qualitative difference between the random effects induced on a wave 
propagating in long range and in rapidly decorrelating random media in time (18 1, for propagation 
media of dimension strictly greater than 1. In fact, it has been shown in 7 that the field c/)^ propagating 
in a rapidly decorrelating medium does not evolve before the scale s = 1 , more precisely the phase and 
the phase space energy evolve at the same propagation scale (s = 1), so that no significant wave 
decoherence can be observed before this propagation scale. In rapidly decorrelating random media 
the scale s = 1 is "universal", in the sense that it does not depend on the statistic of the random 
potential V. 

Theorem |3.1| shows that the phase of the wave field exhibits non trivial stochastic phenomena 
for s — 1/{2kj). As a result, we should expect for larger propagation scale parameters s > 1/{2k^) 
that the random oscillations of the wave field 0^ will evolve faster and faster up to break the wave 
coherence. We show in Section |4] and Section [5] that the loss of coherence appears first on the large 
spatial correlation scales, and then as s increase, it is transmitted to smaller spatial correlation scales 
(see Figure [2]). 
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4 Wave Decoherence for s G (l/(2rt:^), 1) 

This section presents the main result of this paper, it describes the loss of coherence of the wave field 
0e solution of (|7| occurring after the onset of the random phase modulation described in Section [Sj 
On propagation scales e^'^ with s > 1/(2k^) but strictly less than 1, the random phase modulation 
described in the asymptotic e — in the Section [3] begins to oscillate very fast up to break the wave 
coherence, and produce momentum diffusion effect. However, according to |22j the wave decoherence 
does not take place for Sc — 0- To study this diffusion phenomenon we need to consider the scaled 
Wigner transform (19) to capture the wave decoherence for spatial correlation parameters > 0. 



Using the notation introduced in Section pT3| we have the following results. 

In Theorem |4.1[ Theorem |4.2| and Theorem |5.1| we show that the good spatial correlation pa- 
rameter to observe the loss of coherence is 

s, = (l-s)/0, 



so that for Sc < s, we have s > 1/(1 + 0) > 1/(2k^), where 6 is defined by (17) and Kj in Theorem 
3.1 As a result, no decoherence effect can be observed before the propagation scale s = 1/(1 + 6') 
for any spatial correlation parameter Sc 



4.1 



below deals with the case s > 1/(1 + 



< s. Theorem 

for which one can describe the wave decoherence in term of a fractional diffusion, while Theorem |4. 2 1 
deals with the critical case Sc = s = 1/(1 + 9). 



Theorem 4.1 Let us assume that (16) holds. For s G (1/(2k^), 1), and 



< s, 



(28) 



where 9 e (0, 1) is defined by (17), the family of scaled Wigner transform (We)ee{o.i) 



and solution of the transport equation (24), converges in probability on C([0, +oo), (Br, 
to a limit denoted by W . More precisely, for all T > and for all rj > 0, 



defined by ( 19 I 
)) as e 







lim I 



sup dBAW,{t),Wit))>v =0, 
.telo.T] 



where W is the unique solution uniformly bounded in L^(M'^'^) of the fractional diffusion equation 

dtW = -a{9){-^^f'^W, (29) 
P'*) defined by (22). Here, (— Ak)^^^ is the fractional Laplacian 



with 1^(0, x,k) VFo(x,k) e L2(I 
with Hurst index 9 g (0, 1), and 

a{9) = 



2a(0)6ir(l - 9) 



/ d5(u)|ei.u| 



with ei e S''-! and T{z) = /+°° t^^'e-*dt. Moreover, W is given by the following formula 



VF(i,x,k) 



(27r)" 



dqexp(ik.q-a(0)|q|^i)W^o''(x,q), 



(30) 



where Wq stands for the Fourier transform of Wq with respect to the variable k. 

Let us note that we cannot expect a convergence on L^(]R^'') equipped with the strong topology. 
In fact, the following conservation relation |jVFe(t)|j^2(][{2d') = || Wo,e||L2(R2d) is not true anymore for the 
limit W. Moreover, let us note that the Wigner distribution is self- averaging as e goes to 0, that is 
the limit W is not random anymore. This self- aver aging phenomenon of the Wigner distribution has 
already been observed in several studies [2J 31 |SJ [22] for s = 1, and is very useful for applications. The 
proof of Theorem 4.1 is given in Section [7] and is based on an asymptotic analysis using perturbed- 
test-function and martingale techniques. 

Equation ( 29 1 describes the loss of coherence of the field for the particular spatial correlation 



parameter Sc defined by (29) through a momentum diffusion. This fractional diffusion exhibits a 



damping term obeying to a power law with exponent 9 £ (0, 1) describing the decoherence rate of 
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the field wave 4>^. An important point is tliat tlie wave decoherence mechanism is deterministic, it 
does not depend on the particular realization of the random medium. Finally, let us note that W 
does not evolves in x. In fact, in (24 1 the dispersion term k • Vx is of order e""^. As a result, for 
Sc > the dispersion is not captured by our scaled Wigner transform (19) and is small compared 



to the momentum diffusion mechanism. The dispersion effect can be captured by the scaled Wigner 
transform only when 3^ = (see Section [5| . 

The following theorem investigates the special case Sc = s = 1/(1 + ^^), with either 7 > or 
/3 < 1/2. This special case studies the decoherence of the wave envelop 0o itself and not the local loss 
of coherence of the initial condition (see (20 1). The case 7 = and /3 = 1/2, and s = Sc has been 
addressed in Theorem 3.1 since in this particular case 1/(1 + 6) ~ 1/(2ko). Therefore, in this case 



there is no wave decoherence. 



Theorem 4.2 Let us assume that (16 1 holds. For either 7 > or (3 < 1/2, and 

1 



Sr = S = 



1 



where 9 € (0,1) is defined by (17 1, the family of scaled Wigner transform (VFe)eg(o,i) defined by (19 1 
and solution of the transport equation (24l, converges in distribution on C([0, +00), L^(]R^'^)) as e — >■ 
to a limit W defined by 



W is the unique weak solution of the stochastic differential equation 

dW{t,x,k) = - cr{9){-Ai,f^'^W{t,x,k) 
2ia(0) 



(27r)° 



Btidp)e'^ P (ly (<, X, k - |) - M/(t, X, k + |) 



(31) 



with VK(0,x,k) = VFo(x,k) G L'^(R^'^) defined by (23l. Here, {Bt)t is a real Brownian motion on Hg 
the dual space of 

Tie ^ ^ip such that <y3(p) = (/^(-p) and J ^^^_^g |y(p)P < +oo|, 
with covariance function 



dp 

n\d+e 



Moreover, (— Ak)^^^ is the fractional Laplacian with Hurst index 9 G (0, 1), and 

2a(0)6'r(l - 9) 



cj{9) 



with ei e S^'-i and T{z) = /+°° t^-'^e'^dt 



(27r)'* 



/ d5(u)|ei.u|'' 



Let us remark that the proof of the weak uniqueness of (31 1 follows the idea developed in [TBI 
and is the same as the one given in Section |8] to prove that all the converging subsequences of 
(M^c)cg(o,i) have the same distribution. 

This limiting Wigner transform is random because the wave does not propagate enough to observe 



a self- aver aging of the stochastic phenomena. In fact, as shown in Theorem 4.1 for all s > 1/(1 + 6') 



the limiting Wigner transform is self-averaging and is equal to the expectation of the limiting Wigner 
transform obtain in the case Sc — s = 1/(1 -I- 9). 

Let us note that the convergence holds on L^(M^'*) equipped with the strong topology. In fact, 
the conservation relation ||IV£(i)||i2(j|2c!-) = ||M^o,£|!L2(R2d') is preserved for the limiting process W. In 
this scaling the limit is a stochastic process, which is given in the Fourier domain by a random 
phase modulation. As illustrated in Figure [2] the random phase modulation of the Wigner transform 
is caused by the fast phase modulation of the wave field (f)^ itself and described in Theorem |3.1| in 
the asymptotic e -> 0. For s e (1/(2^-^), 1/(1 + 9)) the random phase modulation the wave field (j)^ 
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oscillates very fast up to produce a phase modulation on the Wigner transform at the propagation 
scale with s — + 6). It is the first propagation scale parameter s for which the coherence 
of the wave is affected. Afterwards, for the propagation scale parameter s > 1/(1 + 9) this random 
phase modulation of the Wigner transform oscillates very fast and then average out to obtain the 
determinist wave decoherence mechanism described in Theorem 14. II 

Theorem |4.1| and Theorem |4.2| show a qualitative difference between the random effects induced 



on a wave propagating in long range and in rapidly decorrelating random media in time (see (18l). 
As previously noted, it has been shown in [7] for rapidly decorrelating random media that the field 
4>e does not evolve before the scale s = 1, so that there is no wave decoherence before s = 1. In this 
case the phase and the phase space energy of the wave, which describes the wave decoherence, evolve 
on the same scale. Let us give a probabilistic interpretation to illustrate the difference of the random 
effects caused by the long range and rapidly decorrelating random media. The wave decoherence is 



given by the random operator ( 25 1 , which after homogenization gives rise to an operator of the form 

' J dpaip)[x{k+^)^Xik)). (32) 



Formally, the momentum diffusion of the Wigner is given by the variations of a stochastic jump process 



|34L Chapter 1] with infinitesimal generator given by (32 1, which characterizes its variations. In the 
case of rapidly decaying correlations a{p) E L^(E) so that the variations of the jump process are 
therefore bounded by ©(e^^*) for all s^, that is why we cannot observe wave decoherence phenomena 
in rapidly decorrelating random media for any > 0. The variations of the jump process become 
significant only for the scaling s = 1 and Sc = 0, and that is why we can observe wave decoherence 
only at this scaling. However, for long-range random media the variations of the jump process can be 
large and described by 

J (Mk + ^) - A(k)) e-^-«-a(0) / ^ (A(k + p) - A(k) 

As a result the large variations can balance the small term e^~'^ and give rise to significant momentum 
diffusion. In rapidly decorrelating random media the wave decoherence is only significant for Sc = 
and for a sufficiently large propagation distance, while for the long-range random media the wave 
decoherence occurs first on the large correlation scale and propagates to the smaller ones as the 
propagation scale increase, up to the scaling Sc = and s = 1 (see Figure |2|. We will see in the 
next section that the wave decoherence mechanism in the scaling s — I and ^ for long-range and 
rapidly decorrelating random media is exactly the same, but however, this decoherence mechanism 
has different regularity properties in both cases. 



Let us remark that the spatial frequency of the initial condition (20 1 which is of order e^(''^''<:) is 
low compared to the one of the random medium (~ e"'') on the macroscopic scale e"'*. In rapidly 
deccorelating random media such low frequency sources do not interact with the random medium. 
However, as we have seen in Theorem |4. 1| and Theorem |4.2| low frequency initial conditions interact 
strongly with slowly decorrelating random media. This fact can be useful in passive imaging of a 
target in a slowly decorrelating random medium [19j . 



5 The Radiative Transport Scaling s = 1 and Sc = 

This section describes the evolution of the wave decoherence mechanism in which the momentum 
diffusion and the dispersion are of order 1. The following results have been proved in |22j . but 
we state these results in order to provide a complete and self-contained presentation of the wave 
decoherence mechanism in long-range random media. 

The radiative transport scaling is an important scaling since it is the only one for which wave 
decoherence happens for rapidly deccorelating random media. This scaling provides also wave de- 
coherence for long-range random media, but the momentum diffusion is not exactly the fractional 
diffusion as previously obtained. Even if the decoherence mechanism are the same in long-range and 
rapidly decorrelating random media this result highlights an important qualitative difference between 
the two kinds of random media. 

In the radiative transfer scaling in addition to a momentum diffusion, now we also have a dispersion 
term k • Vx. Using the notation introduced in Section [2.3| we have the following result. 
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(a) 



1/(1 + 0) 






. , / i\ = (1 - s)/e 


Sc = 


/ i \M/a\/e decoherence 




Phase r\iodulatiori\ 


1/(2k^) 1/(1 + 0) i ' 



Se = (1 - s)/0 

Wave decoherence 




1/{2ko) = 1/(1 + 



Figure 2: Schematic representation of the behavior of the wave field 0g solution of Q. s is the 
propagation scale parameter and Sc is the spatial correlation parameter. The phase modulation 
efl^ects appear for s — 1/(2k^) for all Sc < s. Afterward, wave decoherence appears on the low spatial 
correlation scales first, and then propagates to the higher one according to the formula Sc = (1 — s) /9. 
In (a) we represent the behavior of the wave field in the case where 7 > or /3 < 1/2. The dot 
represents the transition between the phase modulation effects and the wave decoherence effects (see 
Theorem 4.2 1. In (b) we represent the behavior of the wave in the case where 7 = and /3 — 1/2. 



Theorem 5.1 Let us assume that (16 I holds. For j > 0, the family {We)ee{Q.i) of Wigner transform 
defined by (19) with s = \ and — 0, solution of the transport equation (24 1, converges in probability 
on C([0, +00), {Br, dg^)) as e to a limit denoted by W . More precisely, for all T > and for all 
rj>0, 

limpf sup dBAW,{t),W{t)) > T]] =0, 
\te[o,T] J 

where W is the unique classical solution uniformly bounded in L^(M^'^) of the radiative transfer equa- 
tion 

dtW + k-\'^W ^ CW, (33) 
with 1^(0, x,k) = W^o(x,k) e ^^(M^'i) defined by ([22|. Here, £ is defined by 



C^{k) = J dpa(p-k)((^(p)-^(k)), 



with LP e C^{W^) and where 



^(P) 



2flo(p) ^ 2a(p) 
(27r)''0(p) (2^)'i|p| 



d\n\d+0 ' 



with 9 G (0, 1). Moreover, W is given by 

I^(i,x,k) = ^ I dydqe^(-^+'^ 'i)e/>-*('>+"'')l?o(y,q + <y), 

where 

vl/(q) = J dpa(p)(e'P 'i-l), 

so that for all to > we have 

VFeC°((0,+oo), fl iJ'X^''')) nL°°([to,+oo), fl i/'=(M2rf)y 



fe>o 



k>0 



(34) 



In Theorem 



5.1 



ffk^^2d--j stands for the fcTh. Sobolev space on E^'', and Wq stands for the Fourier 
transform in both variables x and k. Let us note that the case 7 = 0, has not been addressed in [22], 
because it leads to much more difficult algebra than the cases 76 (0, 1). More precisely, we show in 
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this case the tightness of the family (VFe)^ and show that all the subsequence limits are deterministic 
weak solutions of the same transport equation ( 33 1 with 



£<^(k) 



dpalp — k, 



((^(p) — iy9(k)) with cr(p,(jj) 



2g(p)i?o(p) 

(2^n0(p)2+..2)' 



However, it is difficult to show the weak uniqueness of the limiting transfer equation in the slowly 
decorrelating case. First, the technique used in the proof of Theorem |5 . 1 1 to show the weak uniqueness 
leads to very difficult algebra. Second, it should be possible to use the techniques developed in [TU] , 
but this kind of techniques use a lower and an upper bound of a in terms of |k — p|^(''+®), and we 
just have an upper bound of this form. Nevertheless, we think that the transport equation obtained 
in the case 7 = is still weakly well posed. 

As already discussed in Section 4j we cannot expect a convergence on i^(M^'') equipped with 
the strong topology in Theorem |5.1 since the conservation relation ||VFe(t)||i2(i{2d) = || Wo,e||L2(R2d) 
cannot be satisfied by the limit W. Moreover, the Wigner distribution W^: is also self-averaging as e 
goes to 0, that is the limit W is not random anymore. This self-averaging phenomenon of the Wigner 
distribution has already been observed in several studies [H [H [51 [H] and is very useful for applications. 

Equation ( [33| describes the loss of coherence of the wave field </)<: solution of the random Schrodinger 
equation ([7|. It describes the dispersion phenomenon through the transport term k- Vx, and the wave 
decoherence through the nonlocal transfer operator C defined by ( 34 ) . Finally, the transfer coefficient 
it(p — k) describes the energy transfer between the modes k and p. 

An interesting remark is that the result of Theorem 5.1 does not depend on whether J dpa{p) is 
finite or not. In other words, the radiative transfer equation (33 1 is valid in the two case, slowly and 
rapidly decaying correlations in time (see ([I8])). However, as noted in Theorem 5.1 
in the both cases, behave in different ways. As it has been discussed in Section" pT2l 
rapidly decaying correlations in time, that is / dpa{p) < +00, the radiative transfer equation (33) has 



these equations, 
in the case of 



the same properties as in the mixing case addressed in [3] . In the case of slowly decaying correlations 
in time, that is / dpa{p) — +00, we observe a regularizing effect [221 Theorem 4.1] of the solutions of 
(33 1 which cannot be observed in the case of rapidly decaying correlations in time. This regularizing 



property highlights an important qualitative difference between the cases rapidly and slowly decaying 
correlations, and admits the following probabilistic representation in terms of Levy processes [341 
Chapter 1]. According to ^22j Proposition 5.1], we have 



T4^(t,x,k) = E 



Lgds, k + 



(35) 



where W is the unique classical solution of the radiative transport equation ( 33 1 with initial datum 
Wo G L^(M^'^), and where {Lt)t>o is a Levy process with Lq — and infinitesimal generator (34 1. 



Let us note that (it)t>o is a pure jump process with jump measure a{p)dp. Therefore, because of 
the long-range correlation property in time ( 18 ) the jump process has infinitely many small jumps 
with probability one |34| Theorem 21.3 pp. 136]. This last property of the symmetric process {Lt)t>o 
is the key to show the regularizing effect of the radiative transfer equation (33 1. In fact, according 



to |33L Proposition 1.1], the Levy process {Lt)t>o has a smooth bounded density, which permits to 
obtain smoothness on the variable k. Moreover, the transport term in (33) permits to transfer the 



smoothness from the variable k to the variable x. However, this regularizing effect cannot be observed 
in the case rapidly decaying random media in time. In fact, in this case we have the same probabilistic 



representation (35) but where the Levy process (Lt)t>o has a bounded jump measure cr{p)dp (see 
(I81). In this context, there is a positive probability that on any finite time interval {Lt)t>o does not 
jump |34J Theorem 21.3 pp. 136]. Consequently, in this context there is not enough jump to induce 
regularization on the Wigner transform. 



Let us note that the momentum diffusion in ( 33 1 is not exactly the same diffusion mechanism as the 



one obtained in Section [4] but they are very closed. In fact the momentum diffusion given in Theorem 
4.1| is described in terms of a fractional Laplacian, while in the radiative transfer regime the momentum 
diffusion is described in terms of a nonlocal operator which is not exactly a fractional Laplacian. 
However, this two diffusion mechanisms are anomalous diffusions since they lead to damping terms 
obeying to a power law with exponent 9 G (0,1). We have to wait for a long time of propagation in the 
radiative transfer regime to observe again the momentum diffusion given by a fractional Laplacian. 
This approximation in the radiative transfer scaling is proved in [221 Theorem 5.1] 
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Conclusion 



In this paper we have studied the different behaviors happening on a wave propagating in a random 
media with long-range correlation properties. We have exhibited three different behaviors over a range 
of scales given by Theorem |3.1| Theorem |4.1[ Theorem |4.2[ and Theorem |5.1| These asymptotic be- 
haviors differ strongly from those obtain with the random Schrodinger equation with rapidly decaying 
correlations [3 El E] , for which all the random effects appear on the wave on the same propagation 
scale (s = 1). In the context of long-range correlations, the effects of the randomness appear 
progressively according the scale of propagation. We have seen in Theorem |3.1| that the random 
perturbations induce a phase modulation in term of fractional Brownian motion on the wave field 
itself. This random phase modulation begins to oscillate very fast up to break the wave coherence. 
The first wave decoherence mechanism is described by a random phase modulation on the Wigner 



transform (Theorem 4.2 1 on the large spatial scale. Afterward, this phase modulation on the Wigner 



transform begins also to oscillate very fast up to average out and then gives a deterministic wave 
decoherence mechanism. The wave decoherence first happens on the large spatial scales and then 
propagates to smaller one as the propagation distance increases (Theorem 4.1 Theorem 4.2 Theo- 
rem IKl and see Figure |2]). The wave decoherence mechanism is described in term of an anomalous 



momentum diffusion (Theorem 4.1 and Theorem 5.11 since it obeys to a power law with exponent 
lying in (0, 1). Theorem 4.1 shows that low frequency initial conditions interact strongly with slowly 
decorrelating random media. This result can be useful in passive imaging of a target in a long-range 
random medium |19] . 



6 Proof of Proposition 2.1 



First, to prove ( |13[ ) it suffices to show, for all n > 1 and < ti < • • • < that 

E[F(t„+i, Wih, •),•••, ViU„ •)] = e-3(P)(*-+-*")i/(i„, •), 
on 'H'„i. For n = 1 and < ti < t2, we write 

m,-) -e-0(P)(*^-*i)y(ii,-) + l', 
where Y and V{ti, •) are independent, since they are zero-mean Gaussian variables and 

= J m{dp) </j(p)VKp)(e"B(P)(*^"*^' - e-f(P)(*=^-*i)) 
= 0, 

where tpsip) = e^s^P^''(p(p) and for all G Ti-m x T~Lm. As a result, we have 

E[V^(i2,-)inii,-)] -e-B(p)(*-*^)i/(ti,.). 
Now, let us fix 71 > 2 and assume that for all family (sj)jg{i. such that < si < • • • < s„ 
E[y(s„, Oin^i, •),..., n.s„-i, •)] = e-«(P)(^"-^"-)i>(s„_i, •), 

Then, we write 

y(Wi, •) = e-8(p)(*"+i-*")i/(i„, •) + r, 

where Y and V{tn, ■) are independent as explained above, so that 

E[y(t„+i, Ointi, •),..., V{tr,, ■)] - e-8(P)(*"+i-*")F(t„, •) +E[Y\V{h, •),..., T/(i„_i, •)] 

^ e-8(P)(*"+i-*")F(t„, •) -h E[i>(t„+i, -Wih, •),■•■, l^(<n-i, •)] 
_ e-B(P)(*"+i-*")E[l/(i„, .)|F(ii, •),..., ^(t„_i, .)] 
= e-fl(p)(*"+i-*")i/(t„,.) 

^ (■g-B(p)(tI. + l--^.^-l) _ g-B(p)(^.^ + l-tr^)g-0(p)(^,.-t„-l)^y(•^^^_^^ .^^ 
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which concludes the proof of p3| by induction. Second, to prove (14) it suffices to show, for all n > 1, 
< <i < • • • < tn+i < in+i and {(p, 2p) e Hm X Hm, that 

E[Viin+i){ip)V{tn+i)mV{h, ■),..., y(t„, •)] = E[V{in+i){ip)\Vitu •),•••, V{tn, •)] 

xE[y(t„+i)(V)|T/(ii, •),..., n^n,-)] 

m(dp)</?(p)VKp)(e-s(P)(*"+i-*"+i^ -e-8(P)(*"+i-*")e-s(P)(*"+i-*")). 



For n = 1 and < <i < t2 < t2, we write 

i>(t2, •) = e-B(p)(t2-*i)y(f^, .) + y and i>(i2, •) = e-B(p)(*2-ti)y(i^^ .) _,_ 

where Y and •) are independent as well as Y and V{ti, •), so that 

with 



Now, let us fix n > 2 and assume that for all family (sj)jg{i,...,„} such that < si < • • • < s„ < s„ 
and for all (t/j, ijj) G x 

E[V{sn){^)V{sn)mVisi, •),..., ns„_i, •)] = E[Visn)MV{si, •),•■■, ns„-i, •)] 

xE[F(5„)WI^(si,-),---,^(5n-l,-)] 

+ / m(dp)(^(p)7Kij)(e-8(p)(^""-''") _e-s(P)(^'"-''"-i)e-B(P)(''"-''"-i)). 



Consequently, writing 

V{L+i, •) = e-8(P)(*"+i-*")y(t„, •) + r and V{t,,+,, •) = e-B(P)(*"+i-*") V^(t„, •) + Y, 
where Y and V{tn, •) are independent as well as Y and V{tm •), we have 

E[y(?„+i)(^)y(t„+i)(i^)|1/(ii, •),..., vit„, •)] = E[y(t„+i)(^)|i/(ii, •),..., v{u, •)] 

xE[i>(i„+i)(V')|F(ti,-),---,^(in,-)] 

+ E[fMr(v^)|i/(ii, •),..., n^n-i,-)] 

since 

E[Y\V{h, •),••■, •)] = (e-8(p)(t"+i-*"-i) _ e-s(P)(*"+i-*")e-B(P)(*"-*"-i))y(t„, •) = 0, 

and in the same way E[y|y(ii, •),..., V{tn-i, •)] = 0. Finally, we have 
E[Y{ip)Y{^j)\V{h, •),..., i>(i„_i, •)] = E[V{in+i){ip)V{t^+im\V{h, •),•••, i>(tn-i, •)] 

- E[V{L+l){ip)V{t^){iJt^^,-tJ\V{h, •),.••, ^(tn-i, •)] 

- E[V{t„){^i^^^^jV{t„+,){i^)\V{t,, ■),..., y(i„_i, •)] 

+ E[F(t„)(^,-„^^_,jF(t„)(^i„+,_tJ|i/(ii, •),..., •)], 
where fs{p) — e~3(p)S(^(p) and '0s(p) = e~s*^P'*V'(p)- As a result, we have 

E[Y{^)Y{i;)\V{h, ■),..., i/(t„_i, •)] = Pi + P2, 

where 

A = y(t„-i)(<^f„^,_,„_jv(t„-i)(Vt„+i-t„-J 

- V^(i«-i)(</'t„+i-.t„_jV'(i«-i)(^/'t„+i-t„+t„-t„_i) 

- V'(^«-l)(<^t„+i-i„+t„-t„_JV^(in-l)(V't„+i-t„-i) 

+ V^(i«-i)(<^t„+i-t„+f„-t„_jV^(in-i)(V't„+i-t„+t„-t„-i) 
= 0, 
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and where 



P2 = I m(dp)¥3(p)V'(p)[e 



-B(p)(t>.+l-t.>+l) 



3(p)(tn + l-tTi-l)p-fl(p)(t™ + l-t„-l) 



_ g-0(p)(tr. + l-t„)(g-B(p)(t,. + l-t„) _ g-B(p)(t,x + l-*,>-l)g-fl(p)(tr.-t,>-l)) 



g -a(p)(tr. + l-t,i)(-g-B(p)(tr. + l-t,i) _ g-B(,P)^t™-*,.-ljg-B^P)(t-^ + l-t„-l)) 



-a(p)(tr^ + l-t,^)„-B(p)(^« + l-t>^) 



(1-e- 



-2fl(p)(^„-^„_l)^ 



m(dp)(p(p)V'(p)[e 



-B(p)(t„ 



which concludes the proof of ( 14 ) by induction 



7 Proof of Theorem 4.1 



The proof of Theorem 4.1 is based on the perturbed-test-function approach and follows the techniques 
of the proof of [211 Theorem 2.2]. Using the notion of pseudogenerator, we prove first the tightness 
and then characterize all the subsequence limits. 



7.1 Pseudogenerator 

We recall the techniques developed by Kurtz and Kushner [IHl Section 2.2 pp. 38]. Here, let us 
consider the filtration — Ft/f_ where (J-j) is defined by ( 12 ). Let Ai^ he the set of all J^^-measurable 
functions /(i), adapted to the filtration for which supj^-^E [|/(i)|] < +oo and where T > is 

fixed. The p — lim and the pseudogenerator are defined as follows. Let / and in for all (5 > 0. 

We say that f = p — lim^ if 



supE[|/(t)|] <+oo and lim E[|/(t) - /(t)|] 
t,s *-+o 



Vt > 0. 



We say that / el? (^'^) the domain of A'^ and A'^ f — g ii f and g are in A^*^ and 

'Et[fit+s)]-m 



p — lim 

5^0 



git) 



0, 



where E^ is the conditional expectation given T^. A useful result about the pseudogenerator A"^ is 
given by the following theorem. 

Theorem 7.1 Let f eV{A'). Then 



MUt) = f{t) - /(O) - / A'f{u)du 



is an {Tf) -martingale. 



7.2 Tightness 

In this section we prove the tightness of the family (H/c)(:g(o,i) on the polish space C([0, +oo), {Br, 
according to the following criteria. Following [8, Section 7 pp. 80] and [^T, Theorem 4.9 pp. 62, 
Theorem 4.10 pp. 63], we have the following version of the Arzela-Ascoli theorem for processes with 
values in a complete separable metric space. 

Theorem 7.2 A set B C C([0, +cx)), (6^,(^6^)) has a compact closure ij and only if 

sup dg (5(0), 0) < +00, and \fT > 0, lim sup m7-((7, 77) — 0, 

geB n^Og^B 



with 



mT{g:V)= sup dB^{g{s),g{t)). 

\t-~s\<v 
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From this result, we obtain the classical tightness criterion. 

Theorem 7.3 A family of probability measure (iP'^)^g^Q on C{[0, +00), {Br,dB^)) is tight if and only 
if 



\fr]' > 0, lim^ sup F'{g ; dB,Xg{0),0) > M) = 0, 



!6(0,1) 



and 



yT>0,ri'>0 lim sup F'ig; mT{g,ri) > rj') ^ 0. 

'?^"eG(Oa) 



From the definition (26) of the metric dg^, the tightness criterion becomes the following. 

Theorem 7.4 A family of processes (X'^)jg(Q x) is tight on C {[0, +00), (Br, djs^)) if and only if the 
process ((^% A)^,(„,,))^g(o is tight on C([0,+oo),M) for all X e L^R^'^). 

This last theorem looks like the tightness criteria of Mitoma and Fouque [3TJ[TS]. For any A G X'^(M^''), 
we set 

According to Theorem |7.4| the family {We)^ is tight on C([0, +00), {Br, rfg^)) if and only if the family 
(We,x)e is tight on C([0, +cx)), M), for all A e ^^(M^'^). Furthermore, (T^,), is a family of continuous 
processes. Then, according to [8, Theorem 13.4 pp. 142], it suffices to prove that for all A G L^(M^''), 
the family {Wi^\)e is tight on I?([0, +00), E), which is the set of cad-lag functions with values in E. 
Finally, using that the process takes its values in Br and that the set of all smooth functions with 
compact support in E^'', denoted by C[f (E^''), is dense in L^(E^'*), it is sufficient to show that {We^x)f. 
is tight on X>([0, +00), E) for all A e CJf (E^''). 

Proposition 7.1 For all A e Cg^{R'^'^), the family {W^^x)eeio,i) « tight on V {[0, +00), R). 



Proof (of Proposition 7.1) The proof of Proposition 7.1 consists in applying Theorem 4 pp 
48] which 
let XGCg 



which is based on the perturbed-test-function technique . Throughout the proof Proposition 7.1 
^'^), / be a bounded smooth function, and /o(i) = f{We^\{t)). The pseudogenerator 



A^m) = f'{W,Mt))k^W,Mit) + e^'-'''^^^-'{W.it),CAit)) 



L2(R2d) 



(36) 



where CeX{t,ii,'k) is defined by (25), and 

Ai(x,k) = k- VxA(x,k), 

is well defined since 

supE[||£,A(i)||2 ] < +00. 



Let us remark that (36) blows up as e goes to since s > 1/2. The goal of the perturbed-test- 
function method is to modify the test function /q using small corrector, so that the pseudogenerator 
applied to the test function with its corrector does not blow up anymore (see Lemma 7.2). Let us 
introduce the first corrector 



nW.^xit)) J dxdkiy,(t,x,k) 



+ 00 



(27r)<^i 



/e=gip-x/e° 



X ( A(x, k — 



2e' 



- A(x,k 



P 



du. 



To be able to apply [^Hl Theorem 4 pp. 48] and then prove Proposition 7.1 we have to show the two 
following lemmas. 

Lemma 7.1 VT > 0, and 77 > 

limPf sup \f^{t)\ > 77) = 0, and limsupE [|/i'(t)|] = 0. 

« ^0<t<T ' t>0 

Lemma 7.2 VT > 0, {^'(/q' + fl){t),e € (0, 1), < i < T} is uniformly integrable. 
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The remaining of this section consists in proving Lemma |7.1| and Lemma |7. 2 



Proof (of Lemma 7.1) Using (13 1, we have 



with 



'Cl,eA(t,X,k) 



1 f V{^,dp) 



{2TrYi J 0(p) ~ iefk • p v"^'"'" 2e«<= ^ --y-,-. 

Let us note that the proof of Lemma |7.1| is a direct consequence of the following lemma 
Lemma 7.3 We have 



lime^+'^E 



sup ||/:i.eA(t)|'2 

te[o,T] 



= 0. 



Proof (of Lemma 7.3) First 



ll'^i,eA(i)||i2(R2ci) = J^:^ J ^^cfk 
Let us fixe x, k, and u, and let 



0(p) - ie^k • p 



'(a(x, 



k- — ) - A(x,k 



P 



1 e'P-^A ^ p 



A(x,k- 



2e'* 



According to (17) and since A is a smooth function, we have 0A,x.k,u G ^m- Consequently, 

given 



4>\,x,w,u)-^, ^ is real-valued zero-mean Gaussian process with a pseudo-metric M on [0, T] 



Then, for all (ti,i2) G [0,T]^ 
t 



and 



where 



r7i(dp) 



2(1 - e"B(P)l*i-*2l/'°^"') 

02(p) + e27(k.p)2 



A(x,k-^)-A(x,k+^))^ 



^ r,l*i - *2| [ , a(p) /w , P 
<2^— / dp^— ^ A(x,k- 



2e'* 



A(x,k- 



2e' 



^2j!^^(sup|VkA(x,k)p / dp^^+sup|A(x,k: 



1 



dzamL([0,T]) < 2E[v^2(-|i^)] +2E[y^(-||^) 



2E 

- ^ "^^'P^ fl^(p) + ei(k . p)2 - 2^) - k + ^) 



< 



.(e+2/3)s, 



-62(x,k), 



f 1 /"^/^ 

62(x,k)=/ dpj—^^^ du|VxA(x,k + up 

J|p|<l IPI 



|P|>1 



f)-A(x,k-,|; 



(37) 



(38) 
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Here, diamM{[0, T]) stands for the diameter of [0, T] under the pseudo-metric M. Following the proof 
of [TJ Theorem 1.3.3 pp. 14], and thanks to ([37| and (38 1, we have 



E 



Consequently, 



sup 

te[o,T] 



<Ci 



26(x,k)/ 



^(e + 2/3)sc/2 



/In C2 



T 



dr 



,(s + 7-2/3sc)/2 



ei+^E[ sup \\C^A{t)\\lHR-4 < C / dxdk 
*e[o,T] 



< C'£(s+7-2/3sc)/2 



b(x,k) 



ln(C2^ )dr 



y dxdk6(x,k) j In ^C2-^^dr < +00, 



which concludes the proof of Lemma 7.3 
the proof of Lemma |7.1| □ 

□ 



smce Sr 



(1 - s)/9 and Sc < s < {s + 7)/(2/3), and then 



Proof (of Lemma |7.2[ ) The pseudogenerator A'^(^f(){t) is given by 

+ e'-\f'{w,,^{t)){wMXe{c^A{tW)m^..) 
+ el-^^(w^ea(^))(l^e(^),/:eA(^))^,(„,,)(w^e(^),A,eA(^)) 

SO that one can see that the corrector fl remove the singular terms from A'^ (/q) (t) and we have the 
following result. 



Lemma 7.4 We have 



sup sup e2(i-^)E ||/:,(£i,,A(i))(t)||i2(R2 



< +00, 



and 



sup sup e2(i--)E[||/:^A(t)||2, X ||A,,A(i)||2, ] <+c 
t<Tee(oa) 



Consequently, we obtain 

+ e'-V"(M^.,AW)(W^eW,A,.A(i))^,(„,,^(W^,(i),AA(t))^,(jj,,) 

+ 0(£(l+7)/2+.c)^ 

and supg j E[|^(/q + < +00. That conclude the proof of Lemma 7.2 and then the proof of 

Proposition |7.1[ 



Proof (of Lemma 7.4) A short computation gives 



£,(£i,,A) (i,x,k; 

1 



(27r)2'i 



1 



2e=c 



0(p2)-ie7(k+^) •P2 



(A(x,k- 
(A(x,k 



Pi P2 



A(x,k--P^ + -P^)) 



• P2 ~ - 26"- 26"- 

Pi P2 \ X/ , , Pi , P2 



2e«<= 2e 



-)-A(x,k- 



2e''<= 2e^ 



)), 
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Moreover, the forth order moment of our Gaussian field V is given by 

E[Vitudpi)V{t2,dp2)V*{t3,dp3)V*{U,dpi)] = 

(27r)2^i?(tl - t2,Pl)R{h - t4, P3)'5(Pl + P2)'^(P3 + P4) 

+ {2Tr)^^R{h - t3, pi)R{t2 - U, P3)5{pi - P3)<5(P2 - P4) 

+ {2^f'^R{h - U, Pl)R{t2 - h, P3)(5(pi - P4)<5(P2 - P3), 



so that, using the smoothness of A, (15) and the change of variable p' = p/e'", we obtain 



X C 



'|p|<l \V\ ' ^"'|p|>l 

since s^ — (X — s)/9. In the same way we have 



|^2(R2d) + |lA|1^2(R2d) ) < +00, 



sup sup £2(1 '^)E[||/:eA(t)||i2(R2.) X ||/:i,,A(i)||2,™,J <+oo. 

t<T ee(0,l) 

□ 

□ ■ 

7.3 Identification of all subsequence limits 

In this section, we identify all the converging subsequence limits of the process (VFe)e, given by its 
tightness, as solutions of a deterministic diffusion equation. For the sake of simplicity we still denote 
by (We)e such a subsequence. Let us note that in this case all the limit processes are therefore 
deterministic, and the convergence of the process (T^e)e also holds in probability. We will see that 
this fractional diffusion equation is well posed. In particular, this will imply the convergence of the 
process (VFe)e itself to the unique solution of this diffusion equation. 

Proposition 7.2 Let W be an accumulation point of {Wt)^. Then, W is the unique strong solution 
of the diffusion equation 

dtW ^ ~a{e){-^^f'^W, (39) 



with Ty(0,x,k) = Wo(x,k) defined by (22), and 



iT{e) = — / dS'(u)|ei • u| . 

To prove this proposition we use the notion of pseudogenerator introduced in Section |7.1| and 
the perturbed-test-function technique that we have already used in Section |7.2| for the proof of the 
tightness. Thanks to the pseudogenerator we are able to characterize the accumulation points of (We)e 
as solutions of a well-posed martingale problem, but as we will see it will turn out that the martingale 
problem is only a deterministic partial differential equation. However, as we saw in Section |7.2| the 
pseudogenerator yl'^(/o) has singular terms, so that we have to modify the test-function /g with a 
small corrector ff in order to remove these singular terms. As a result, A'^{fQ + fi) has no more 
singular term but it is not yet a generator of a martingale problem. Consequently, we will introduce 
a second small corrector /| so that A'^{fQ + /i + /I) is now a generator allowing us to characterize 
the accumulation points. 



Proof (of Proposition 7.2) In this proof we use the following notation 

(fi (g) ■0(xi,ki,X2,k2) = (^(Xi,ki)-0(x2,k2). 

Let us introduce the second corrector 

m) -el-V'(M^.,A(i))(W^eW, i/l,e(t))^2(R2.) 

' nW,^x{t)){W, it) ® it) ,H2At)) ^2 , 
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where 



i?i,,(t,x,k) 



+00 



du 



X gi(Pl+P2)-x/c°gj(ti-t)(pi+P2)-k/e° 



B(p2)-ie7(k- ^) - PS 
1 



II (eI [V{u/e'+\ dp,)V{u/e'+\ dp^)] -E[V{0, dpi)V{0, dps)] 



X, k 



Pi P2 



A(x,k 



Pi , P2 



0(P2) -ze^(k+ ^) - Pa 



2eSc 2e«-' ^ 2e''- 2e 

A(x, k + - P^) - A(x,k + + P^) 



and 



i/2,e(t,Xl,ki,X2,k2) = 

1^^ du II {E\[V{ule'+\dp,)V{u/e'+\dp2)]~E[V{Q,dp^)V{Q,dp2)\ 



(2t 



X gipi-xi/c''gip2-X2/e''gi(u-t)(pi-ki+p2-k2)/c° 



1 



0(Pi) - ie'^ki • pi 



X (A(x2,k2 - ^) - A(x2,k2 



2e''<= 

P2 

2e'*<= 



Let 



r(i) = /oW + /rw + /2W- 



According to Theorem 7.1 (Af^, (t))^^^ is an (J^j )-martingale. That is, for every bounded continuous 
function <&, every sequence < si < • • • < s„ < s < t, and every family (/ij)je{i....,n} G L^^jjSd-jn^ 



have 



E 



0, 



(40) 



where after a long but straightforward computation using ( |14[ ), we have 

^'(/o + II + /2)W =e^-V'(M^ea(t))(l^.(t), Gi,,A)^,(„,,^ 

+ e^-'f"{W,^x{t)){W,{t) ^ T^.W, G2,eA)^,(^,,) 



(41) 



0(1), 



with 



Gi,,A(x,k) = 



(27r)° 



dpi?o(p) 



— —7^ ^-^^ (a(x, k) - A(x, k - 

^ ^(x,k+^)-A(x,k) 



and 



G2,eA(xi,ki,X2,k2 



0(p)-ie7(k+2fb) -P 

g{p)Ro{p)e'P<^^-^^y' 



(42) 



dp 



{2n)d J ^ (g(p) - ie7ki • p)(2g(p) - ie^{ki - k2) • p) 



A(xi,ki 



2e^ 



A(xi,ki 



P 



(43) 



A(x2,k2-— )-A(x2,k2+ — 



Let us remark that the second corrector /| is effectively small. 
Lemma 7.5 

limsupE[|/|(t)|] =0. 



23 



As a result, using (41 1, Lemma [7. 1[ Lemma 7.5 we have the following result 
Lemma 7.6 

/■* 

M 



is a martingale where 



GiA(x,k) = -fT(0)(-A)«/2^(x,k). 



More particularly, let us consider f be a smooth function so that f(v) = v, for all v such that 
\v\ < r||A||i2(R2d), then 



is a martingale with a quadratic variation equal to 0. Consequently, Mx = 0, that is W is a deter- 
ministic weak solution of the diffusion equation (39 1. 

To show the weak uniqueness of this equation, let us assume that Wq — 0. Moreover, it is clear that 
this diffusion equation with initial condition Aq G L^(M'*) admits a unique strong solution that we 
denote by X" of the form (ISOl). As result, for all T > 



{W{T),Xo)^.,^.,.^{W{T),X'{T)) 



{W{t),dtX\t)+ai0)i-Af^'X'>it)) 



dt^O, 



for X^{t) — X^{T~t). As a result, all the accumulation points are also strong solutions of the diffusion 



equation (39 1. That concludes the proof a Proposition 7.2 



Proof (of Lemma 7.5) According to (14 1 
t 



V{u+^,dp,)v{u+-^,dp: 



gS+7 

= e 



E 



ViO,dpi)V{0,dp2) 



Consequently, 



mt) = f'{W,Mt)){WM,Hi,,{t))^ + f{W,,x{t)){W,{t) W,{t),H2,e{t)) 



L2(R4d) 



where 

i7i,,(i,x,k) 



(27r)2'ii2 



y(t/e^+T,dpi)l/(t/e^+T,dp2 



o«(Pl+P2)-x/£ = 



(fl(Pi) + fl(P2) - ie'^k • (pi + P2)) 



fl(p2)-ze^(k-^).p2 



0(p2)-ieT(k+^) ■p2 



A(x,k 

A(x,k 



Pi P2 



2e^<= 2e^' 

Pi P2^ 

2e^<= 2e''<= 



XU,k- 



Pl , P2 

2e' 



XU.k- 



- + —)) 

Pi J>2^ 



(27r)" 



1 



20(P) L0(p)+*e7(k- Jb)-P 



A(x,k--^)-A(x,k) 



g(p) + ie7 k + . p 



A(x,k)-A(x,k+^) 
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and 



-ff2,£(^,Xi,ki,X2,k2) 



'(27r 



i2d,-2 



V{t/e'+'^,dp,)V{t/e'+'^,dp2) 

pipi-xi/c" ip2-X2/e'' 



(0(Pi) - ie^^i ■ Pi)(0(Pi) + 0(P2) - ie''(ki • pi + k2 • P2)) 
^ (Mxi,ki-|^)-A(xi,ki + |^ 

X (A(x2,k2-U:)-A(x2,k2 + |^ 
~j2nrJ P (0(p) - ^e^ki . p)(20(p) - le^k, - k^) ■ p) 

P \ ^^ , , P 



A(x2,k2 



2e'* 



In the same way as in Lemma |7.4[ we have 

SUpE[||iJi,,(i)||2 +||iJ2,,(t)|| 



A X2,k2 



< g-2(e+2/3)se 



26" 



X C 



X (l|-C'xA||^2(Jl2d) + II VxA||^2(R2d) + ||A||^2(R2d) j , 



(44) 



and 



1 + 7 - (0 + 2/3)s, > ^ s > (1 - 70/(2/3))/(l + 0/(2/3)) = 1/(2k^), 
that concludes the proof of Lemma |7.5[ □ 



Proof (of Lemma 7.6) We want to pass to the hmit e ^ in (40 1. Making the change of variable 
p' = p/e*'= in (42 1 and (43), the e-power in (41 1 becomes 1 — s — 9sc- Consequently, the particular 
choice Sc — {1 — s)/6 is made to obtain a nontrivial limit of (41 1 as e —> 0. Consequently, with the 
change of variable p' = p/e^" in (43 1, and using the Riemann-Lebesgue Lemma and the dominated 
convergence theorem, we have 

lim||ei-"G2,eA||22(g4d) =0, 



since Sc — {I ~ s) /9 < s. Moreover, with the change of variable p' — p/e*" in (42l and the dominated 
convergence theorem, we have 



lim||ei-^Gi,,A-GiA||i2(R2.) = 



(45) 



where 



GiA(x,k) 



dp- 



\d+0 



A(x, k + p) + A(x, k - p) - 2A(x, k) 



«(0) 

= -(7(0)(-A)^/2^. 

Consequently, we have 

A'{r^ + n + = r{w,^x{t)){w.xt), gia)^2(r2.) 

that concludes the proof of Lemma |7.6| using Lemma |7.1[ Lemma |7.5[ □ 



(46) 



0(1), 
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8 Proof of Theorem 14.2 



The proof of this theorem is quite similar to the previous one. The proof of the tightness is exactly 
the same. However, in this theorem we have assumed Sc = s, but also either /3 < 1/2 or 7 > 0. To 
characterize the accumulation points we use exactly the same perturbed test functions as in the proof 
of Theorem 14.11 to obtain 



(47) 



0(1), 



where Gi^^X is defined by (42), and G2,eA by (43) with s — Sc- However, we still have (45) for the 
drift term /', but with the change of variable p' — p/e* there are no fast phases anymore, and we 
obtain that 



Mf^x{t) = f{Wx{t))~f{Wxm~ / f'{Wx{u)){W{u),Gi\)^,^^,,^ 



+ f"{Wx{u)){W{u) ® W{u),G2{X, A))^,(„,,)du 



is a martingale where Gi is defined by (46), and 



G2(Ai, A2)(xi,ki,X2,k2) 



1 

{2nY 



dp 



,ip-(xi-X2) 



X (Ai(xi,ki-|)-Ai(xi,ki + |; 

A2(x2,k2 - ^2(x2,k2 + ^] 



(48) 



In this theorem the second order term /" has not been killed by the Riemann-Lebesgue Lemma, 
so that the limiting point W is not deterministic anymore. As a result we need to study the finite 
dimensional distributions 

limE[n(W^e(t,),A,);^,(„,,; 

to characterize all the accumulation points of (W^)^. To do that, we follow the technique used in 
[iniES] and let us consider the tensorial process W^^ (t) = (g)^! W<:(t) on L2(R2dM)^ ^^^^ 
way as the case Af = 1, we can show that W^'^ {t) is a tight process in L^(M^'^*^) equipped with the 
weak topology and for all its accumulation points VF*^ 



Ahit) = W^\t) - W-f (0) - / (iy*'^(u), Gf A + Gf A)^,(j,,,,,)d^i 



is a martingale for all X E L 



2(m2dM 



and 

G2'-^A(xi,ki, • • • ,XM,kM 

M 

E 



^ {2t:Y 



). Here, Gf and G2 are defined by 

M 

G','X = -J2a{e){-A^//'X 



dp 



«P-(Xjl-Xj2) 



X (^A(xi,ki,--- ,Xji,kj, - ,XA/,kA/) - A(xi,ki,--- ,Xji,kji 

X (^A(xi,ki,--- ,Xj2,kj2 - ,XM,kM) - A(xi,ki,- 



! ^32 , kj2 



P 

2' 
P 

2' 



,xjv/,kM, 

,XM,kM, 
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As a result, E[Vl^ ] is a weak solution of the differential equation 



dt\^' = (Gf + )A*^ 



(49) 



Let Aq G i^(E^''*^) such that its Fourier transform with respect to (ki,--- ,kA/), Aq, belongs to 



C^(^2dMy Solving (|49|) in the Fourier domain, we show the existence and uniqueness of a smooth 



function A of (49 1 in the strong sense with initial condition Ao- As result, if E[iy (0)] = 0, for all 



T > 0, 



L2(R2dM) 



dt = 0, 



for \^\t) = \^'\T — t). Consequently, by a density argument, we obtain the weak uniqueness of (49l 
and the moments 



E 



W*^(t,xi,ki,... ,XM,k 



M) 



are therefore uniquely determined for all accumulation point W. Let 



1 



W{t,^,k) ^ j—j^ I dqW^o''(x,q)exp(zk.q + z / St(dp)e'P-^(e 



Using the ito's formula and the weak uniqueness of (49 1, we obtain for all accumulation point W 
of {W,^h 



M 



E 



M 

M^*'^(i,xi,ki,--- ,XM,kA/)] =E[niy(i,x„k,) 



Consequently, we have identified the one-dimensional finite distributions for all accumulation point 



N N N 

limE[n (W-.W, A,);^.(„.,J = E[n (W^fe), A,)2^.(„.,J = e[ J] (^t), A,);^.^^,.,) 



3 = 1 



To conclude the proof of Theorem [42] following a classical argument regarding the proof of unique- 
ness of martingale problems [211 Proposition 4.27 pp. 326]: If the one-dimensional distributions of 
two solutions are the same, then their finite dimensional distributions are also the same. Consequently, 

N N N 

limE[n (M^.fe), A,)2^.(„.,J = E[n A,)2^.(„.,J = E [ J] (W^fe), A,)2^.(„. 
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